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Many stars are surrounded by disks of dusty debris formed in the collisions of 
asteroids, comets and dwarf planets. But is gas also released in such events? 
Observations at sub-mm wavelengths of the archetypal debris disk around (3 
Pictoris show that 0.3% of a Moon mass of carbon monoxide orbits in its debris belt. 
The gas distribution is highly asymmetric, with 30% found in a single clump 85 AU 
from the star, in a plane closely aligned with the orbit of the inner planet, (3 Pic b. 
This gas clump delineates a region of enhanced collisions, either from a mean 
motion resonance with an unseen giant planet, or from the remnants of a collision 
of Mars-mass planets. 


significance level of 4a, suggesting that 
it may influence the innermost dust 
distribution. 

ALMA also detected the disk in the 
12 CO J = 3-2 transition (Fig. IB), with a 
clear velocity gradient along the major 
axis, illustrated in the position-velocity 
(PV) diagram (Fig. 2). This shows the 
characteristic distribution of a broad 
belt of orbiting gas, with inner and 
outer radii of 50 and 160 AU, and a 
peak around 85 AU. No gas emission is 
seen inside 50 AU. The CO distribution 
is, on average, a factor of 2 brighter to 
the southwest, and a similar asymmetry 
was seen in the mid-infrared emission 
(72). Unlike the sub-mm continuum, 
the CO clump is offset by -5 AU above 
the main disk plane (Fig. IB) and is 
more closely aligned with the inner 
planet b Pic b and a secondary disk 
seen in scattered light (5, 9). 

The CO J = 3-2 emission totals 7.6 
± 0.8 x 10- 20 W/m 2 . Assuming a range 
of possible excitation temperatures 
from 20 K [measured from UV absorp- 
tion lines (13)], to 85 K [the dust equi- 
librium temperature (10)], this 
corresponds to a CO mass of 1.7 x 10 20 
kg (0.0023M Moon ) within a factor of 2. 
Toward the star itself, the total CO 


column is 2.5, 


(+2.5, -1.2) 


10 1 


Debris disks are the end product of collisional cascades of km- sized 
bodies orbiting stars (including the Sun), and are normally thought to 
contain negligible gas. At a distance of 19.44 parsec (7) and age of 20 
million years (2), (3 Pictoris is one of the closest, brightest and youngest 
examples. Its edge-on disk was the first to be imaged in scattered light, 
showing the distribution of micron-sized dust grains (3). Various subse- 
quent observations have provided evidence of infalling comets within a 
few Astronomical Units (AU) of the star (4), a massive planet at -10 AU 
(5), as well as atomic gas extending out to -300 AU (6); it is still unclear 
how these features are linked. By observing such debris disks at mm 
wavelengths, it is possible to trace the mm-sized dust and, by inference, 
the parent bodies of the collisional cascade, known as planetesimals (7). 

We observed (3 Pic using the Atacama Large Millime- 
ter/submillimeter Array (ALMA) at a wavelength of 870 pm in both the 
continuum and the J = 3-2 12 CO line with a resolution of 12AU (8). The 
continuum image (Fig. 1 A) shows maxima in the surface brightness -60 
AU either side of the star. At separations from 30-80 AU the continuum 
is, on average, 15% brighter to the southwest. These results indicate that 
the mm grains he in a broad slightly asymmetric belt nearly co -located 
with the disk of sub-micron reflecting dust (9). The total flux of 60 ± 6 
mJy corresponds to a dust mass of 4.7 ± 0.5 x 10 23 kg (6.4M Moon ), if we 
assume a standard dust mass opacity (0.15 m 2 kg -1 at 850 pm) and tem- 
perature of 85 K (10). 

The planet p Pic b was close to maximum SW elongation at the time 
of the observations, with a projected separation of 8.7 AU (77). Its loca- 
tion is coincident with a dip in the continuum surface brightness at a 


suits from UV absorption line meas- 
urements range between 0.6-2. 1 x 10 15 
cm -2 (13), indicating that at least 10% 
of the CO along this line-of-sight lies in 
front of the star. 

The attenuation through the disk due to dust or CO self-shielding is 
low and, at radii >50 AU, CO is destroyed mostly by UV photons from 
the ambient interstellar medium (ISM) (14). The CO photodissociation 
timescale in the unshielded outer disk will be -120 yrs (75) - substantial- 
ly less than the 600-year orbital period at 85 AU. Unless we are observ- 
ing p Pic at a very unusual time (i.e., <100 years after a collision) then 
the CO must be continuously replenished, with a steady-state production 
rate of -1.4 x 10 18 kg/yr. The source of CO is likely the icy debris 
(grains, comets and planetesimals) in the disk (16). In our own solar 
system, comets are mostly composed of refractory silicate grains togeth- 
er with H 2 0 and CO ices, with a typical CO/H 2 0 ratio of 0.01-0.1 (77). 
Although CO sublimates at -20K - substantially less than the 85 K dust 
temperature in the P Pic disk - it can be trapped in H 2 0 ice at tempera- 
tures up to -140 K (18), which could occur in radiative equilibrium as 
close as 30 AU from p Pic (14). At larger distances, photodesorption or 
collisions can release this CO into the gas phase. Photodesorption (a.k.a. 
UV sputtering) of mm-sized ice grains (19) in a dust clump with 10% of 
the total disk mass would produce CO at a rate that is -10 times lower 
than our observations show. The CO linewidth measured directly toward 
the star is -1 km/s (Fig. 2), suggesting the velocity dispersion is lower 
than the 6km/s thought necessary for collisions to directly vaporize ice 
(20). However icy parent bodies being shattered in multiple lower veloc- 
ity collisions may release much of the entrapped CO as gas, leaving the 
less volatile H 2 0 ice behind. To sustain the gas production rate with a 
CO ice release fraction of 0. 1 , the required mass loss of colliding grains 
would be -1.4 x 10 19 kg/yr, equivalent to the destruction of a large com- 
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et every 5 min. If sustained over the lifetime of P Pic, ~50M Earth of icy 
bodies will have been removed. 

Continual replenishment of CO may explain another puzzle. In the 
absence of a braking mechanism many of the atomic species such as Na 
found around p Pic should be blown out from the central star under the 
force of radiation pressure (21-23). Carbon is overabundant relative to 
other metals, by factors of 18 to 400 (27, 22), and has been proposed as 
the mechanism that allows the ionized gas to remain bound to the star. 
C+ couples to other ions through Coulomb forces; and since it feels 
weak radiation pressure from the central star, it tends to brake all the gas. 
The origin of the unusual carbon overabundance may be explained if CO 
is being rapidly released in the clump, then dissociated and ionized be- 
fore being spread throughout the disk. In steady state, the relative C+/CO 
abundance implies that C+ is being removed on a timescale of 3 x 10 3 - 
10 4 yrs; since this is several orbital timescales at 85 AU, the braking gas 
can be spread throughout the disk. 

If the gas is on circular Keplerian orbits, each point in the PV dia- 
gram corresponds to one of two points in the disk. Figure 3 shows possi- 
ble deprojections of the data to form the face-on gas distribution (8). A 
compact clump of CO is found at 85 AU radius which contains -30% of 
the flux, along with an extended “tail” of emission. The radial distribu- 
tion of the mm dust can also be derived (8, 12), and shows excess around 
this radius in the southwest (Fig. 3C). The mm dust shows an abrupt 
drop in density around 130 AU radius, coincident with a change in slope 
in the radial distribution of scattered light. This was interpreted as a tran- 
sition from a dust-producing planetesimal belt to a region dominated by 
submicron- sized grains blown out by radiation pressure (9, 24). 

What is the origin of the CO clumps and corresponding mid-infrared 
and continuum features? There are two possible interpretations, both 
requiring planet-mass objects in the outer disk, but with different predic- 
tions for the deprojected gas distribution. In the first (Fig. 3 A), outward 
migration of a planet traps planetesimals into both the 2:1 and 3:2 mean 
motion resonances, resulting in a distribution with two clumps of asym- 
metric brightness on opposite sides of the star that orbit with the planet 
(25). Planetesimal collisions would occur most frequently in the clumps 
which would then be the most vigorous production sites of both CO and 
the short-lived micron-sized grains seen in mid-infrared images (72). 
The morphology of the clumps constrains the planet’s mass and migra- 
tion parameters (25). Notably the planet must be >10M Earth to have cap- 
tured material into its 2:1 resonance, and it would currently be close to 
the inner edge of the gas/dust belt, -90 degrees in front of the SW 
clump. The CO clump’s “tail” would point in the opposite direction to 
the rotation, and its’ length given by the ratio of CO photodissociation 
lifetime to the synodic period. 

The alternative interpretation (Fig. 3B) is that the CO and micron- 
sized dust originate in a single recent collision (72, 26). Since the clump 
also contributes -10% of the flux in our sub-mm continuum image, the 
parent body must have -Mars mass, as giant impacts typically release 
-10% of the progenitor's mass as debris (27). While collisional debris 
persists in a clump for -100 yr after such an event, it is more likely that 
the collision occurred -0.5 Myr ago, with the asymmetry arising because 
the orbits of all collisional debris pass through the collision point, a re- 
gion which would retain a high density and collision rate (26). Such a 
clump would be stationary, and the CO “tail” would lie in the rotation 
direction, with a length given by the ratio of photodissociation time to 
sidereal period (Fig. 3B). 

Tentative evidence has been published for the clump’s motion (28), 
which if confirmed would favor the resonance interpretation. Either way 
these observations provide a valuable opportunity to ascertain the preva- 
lence of planet-sized objects at large orbital distances in this disk. 

The CO gas detected in p Pic is unusual among debris disks. So far, 
two others have been detected in CO [49 Cet and HD21997 (75)] and 
one in OI [HD 1725 5 5 (29)]. Most debris disks are presumed to contain 


icy grains, so we would expect all such systems to host CO and its pho- 
todissociation products. The high collision rates found in clumps may 
substantially enhance the gas abundance. If CO is formed in collisions, 
its brightness will scale as the collision rate, which is proportional to v r 
n d 2 , where v r is the collision velocity and n d is the number of dust grains. 
Both n d and v r are enhanced in clumps due to resonances and giant im- 
pacts, increasing collision rates by factors of up to 10-100 (7, 26). The 
presence of these clumps can increase the overall collision rate and CO 
flux from a disk by an order of magnitude. The CO and compact clump 
in the p Pic disk indicates a period of intense activity, driven by planets 
or planet collisions. 
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Fig. 1. Images of p Pic using ALMA in continuum and line emission. Both 
have been rotated by +29 degrees, the position angle of the main dust disk. 
The beam size (0.56 x 0.71 arcsec) is shown lower left, and the locations of the 
star and planet on the date of observations are indicated by the asterisk and 
the cross. (A) Continuum emission at 870 pm from the mm dust; the rms noise 
level is 61 pJy beam -1 . (B) Total J = 3-2 CO line emission (rest frequency 
345.796 GHz); the noise level is 0.02 x 10“ 20 W rrf 2 . The planes of the main 
and secondary dust disks (9) are shown by the dashed lines. 
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Fig. 2. Velocity distribution along disk major axis. This 
shows the observed velocity of the CO emission relative to 
the velocity measured toward the star (which has a centroid 
of +20.3 ± 0.1 km/s, in the barycentric reference frame). Gas 
on the brighter SW side is approaching us; the two dashed 
white lines represent the velocity we would measure from 
gas at the inner and outer radii of the CO belt. The solid red 
lines show the true orbital velocity, the maximum observed 
velocity from gas at the tangential point, as a function of 
distance from the star, assuming a stellar mass of 1 .75M SO iar- 
Contours are 10, 20, 30, 40, 60, 80% of the peak (0.06 
Jy/beam). The spectral/spatial resolution is illustrated lower 
left. 
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Fig. 3. Deprojected distributions of gas and dust around p Pic (8). (A and B) Two 

possible face-on distributions of the CO assuming circular orbital motion. (A) The two- 
clump distribution, interpreted as mean-motion resonances with an inner planet. In this 
case the bright southwest resonance lies on the far side of the star and is approaching 
us. (B) Single clump distribution, interpreted as a collision of ~Mars-mass objects; in 
this case the southwest clump is stationary and the “tail” of CO points in the direction of 
rotation. (C) Radial distribution of mm-sized dust in the southwest and northeast sectors 
of the disk, obtained by fitting consecutive annuli to the continuum distribution on the 
two sides (12). The clashed line is an axisymmetric model for the parent body 
distribution predicted from scattered light images of the sub-micron grains (24). Surface 
density units are cross-section area per unit disk area, in AU 2 /AU 2 . 
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